Spin- Charge Locking and Tunneling into a Helical Metal 
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We derive a kinetic equation for the electrons moving on the surface of a three-dimensional 
topological insulator. Due to the helical nature of the excitations backward scattering is suppressed 
in the collision integral, and the spin dynamics is entirely constrained by that of the charge. We 
further analyze the tunneling between the helical and a conventional metal or ferromagnet. We 
find that the tunnel resistance strongly depends on the angle between the magnetization in the 
ferromagnet and the current in the helical metal. A nonmagnetic layer on top of the helical metal 
amplifies the current-induced spin polarization. 



Topological insulator*^ have recently attracted con- 
siderable interest, especially after their experimental dis- 
covery in twcJS and three dimensioniPH3 While insu- 
lating in the bulk, such materials possess gapless helical 
edge states whose existence depends on - and is protected 
by - time reversal imparlanc e 8 ™ 14 !. This makes the lat- 
ter robust against time-reversal symmetric perturbations 
(such as impurity scattering) and at the same time very 
sensitive to time-reversal breaking ones (such as mag- 
netic fields). When the topological insulator is a three- 
dimensional system, the gapless excitations are confined 
to its surface and form a two-dimensional conductor and 
presents novel and interesting properties, see for example 
Ref. |15| for a recent summary. In particular, Burkov and 
Hawthorn^ considered the problem of spin-charge cou- 
pled transport on a helical metal and derived diffusion 
equations for charge and spin. They predicted a distinc- 
tive magnetoresistance effect when the helical metal is 
placed between a ferromagnet and a normal metal. In 
this paper we extend their work in several ways. We first 
derive a kinetic equation which is valid even beyond the 
diffusive regime. In this latter regime we obtain a dif- 
fusion equation which agrees with that of Burkov and 
Hawthorn as far as the charge component is concerned. 
On the other hand for the spin density we find a differ- 
ent behaviour, namely the spin dynamics is constrained 
to follow the charge one. Secondly, we consider the effect 
of bringing the helical metal in contact with a ferromag- 
net and discuss its unconventional magnetoresistance. 

For the simplest case the effective Hamiltonian describ- 
ing the s urface states of a topological insulator has the 
forrrPia 



H = i>j?k x e 2 ■ cr, 



(1) 



where the parameter vp is the velocity of the gapless ex- 
citations, e z is a unit vector perpendicular to the surface, 
k is the two-dimensional momentum operator, cr are the 
Pauli matrices, and units of measure such that h = 1 
have been used. The eigenstates of H form two bands 
with linear dispersion, e± — ±i>i?fc, and we will assume 
in the following that the Fermi energy is located deep 
enough in the upper band for states in the lower band 



to remain fully occupied and thus not relevant for the 
dynamics of the system. 

Because of the helical nature of the excitations, such a 
surface conductor is called a helical metal and presents 
novel and interesting properties. For instance the veloc- 
ity operator is given by x = vf^z x c, so that the particle 
current becomes j = 2vf&z x s, s being the spin polar- 
ization. This means that the particle current is entirely 
constrained by the spin density or that, vice versa, the 
in-plane components of the spin density are constrained 
by the particle current. 

Such a constraint will also become apparent later in 
the kinetic equation for a disordered helical metal. The 
equation will be valid when the Fermi energy is far from 
the Dirac point and satisfies the condition £f » 1/r, 
with r the scattering time. We will follow a procedure 
similar to what was done for graphene in Ref. [18] . The 
starting point is the retarded Green function which in 
the absence of disorder reads 



where 



and 



with 



G SS ' = G [a ] ss > + G • [<t] ss i, 



G« = ± (G% + G«) 



G^^kx e z (G» - G R ) , 



G± = (eTv F k + ^ + 10+) 1 



(2) 
(3) 

(4) 
(5) 



and k being the unit vector in the k-direction. For clarity 
we included in Eq. ^ the spin indices s and s'. The 
Green function G has the two-by-two matrix structure of 
the Keldysh formalism 



G(R,T;k,e) 



G R {R,T;k,e) G(R,T;k,e) 
G A (R, T; k, e) 



(6) 



where R and T are the center-of-mass space and time 
coordinates, while k and e are the Fourier transformed 



2 



variables of the relative coordinates. The left-right sub- 
tracted Dyson equation reads 

d T G +^ {e z x a ■ d R ,G} +i[H,G] = -i[t,G], (7) 

where [,] and {, } are the commutator and anticommu- 
tator. On the right-hand-side of the equation the self- 
energy E appears. For a delta-correlated impurity po- 
tential with (V(x)V(x')) = M 2 <5(x — x') and within the 
Born approximation the self-energy is proportional to the 
Green function, integrated over the momentum 



E(R,T;e) 



d 2 k 
(2tt) 2 



G(R,T;k,e). 



We define a quasiclassical Green function as 



ff(R,T;k,e) = - / d£G(R, T; k, e), 



(8) 



(9) 



where £ = vpk — fi and the integration is performed in 
the vicinity of the Fermi energy. From Eqs. |2|-([5|, one 
obtains the retarded component of g for energies close to 
the Fermi level (|e| <C vpkp) as 



1 1 



9 = 



k x e z ■ cr, 



(10) 



i.e. g R is a projector on the upper band of the Hamilto- 
nian 0. The self-energy reads 



£ = — (9), 

T 



(11) 



where 1/r = nN u 2 , N — kp/(2irvp) is the single- 
particle density of states at the Fermi energy and (g) 
is the average of the quasiclassical Green function over 
the Fermi surface. Finally, for the Keldysh component 
of the Green function we find from Eq. ^ the kinetic 
equation 

drg+^-{e z x cr-d R ,g} 



ivpk} 



kxe, 



1 1/ i 

-g+ -(g) 



^{kxe 2 -(T,( 9 )} 



(12) 



T 27 

Notice that this is a matrix equation in spin space, 
9ss' = <?o[co]ss' + g • [o"]ss'- However, its structure can 



be considerably simplified. Equation (12 1 is derived un- 



der the assumption that vpkp is much larger than all 
other energy scales in the problem. Therefore, the lead- 
ing contribution to g must commute with k x e z • cr, or in 
other words, g is diagonal in the eigenstates of H. Since 
the lower band has no density of states at the Fermi level, 
only the upper band contributes to g which is then pro- 
portional to the upper-band projector, i.e. 



9ss' = go ( Wo]ss> + k x e z ■ [cr]. 



(13) 



From the spin-trace of Eq. ( 12 ) we obtain then a kinetic 
equation for go, 

— W(tp - <p') [g o (cp) - go (</)] , 

(14) 



which is just the standard kinetic equation for a metal 
with an angle dependent scattering potentia l 19 * 20 !. Here, 
starting from a short range potential, we find 



1 



1 



WUp, <p') = -(1 + k • k') = -(1 + cos(<£ - if/)), (15) 

T T 

where the cosine term accounts for the absence of 
backscattering. The spin-dependent contributions to g 
can be reconstructed from go, cf. Eq.(13|: 



g x « k y g , and g y m ~k x g . 



(16) 



The g z component is nonzero only to subleading order in 
1/vpkp and is after some algebra determined as 



vpki 



-r(ky9o) 



k'l, 



(k x go) + v F k y d x go - vpk x d y g Q 
(17) 



Equations ( 14 ), ( 16 1 and ( 17 1 are one of the main results 



obtained in this paper. 

It is instructive to study the angular average of the 
kinetic equation, since the latter is closely related to the 
continuity equation for the observables. Recall, for ex- 
ample, that the particle and spin densities are obtained 
from the quasiclassical Green function by taking the an- 
gle average and integrating over the energy^, 



No 
2 



de(.go) + N e(/), 



No 
' 4 



Ms), (is) 



where <j> is the scalar electrical potential. From Eq. (14) 
we obtain the continuity equation for the density as 



drigo) + d R ■ (vpkg ) = 0. 



(19) 



Using Eq. (16 1 we find Vp(k x go) — —vp(g y } and 



vp(k y go) — vp(g x ), i.e. we verify the general relations 
between the particle current and the in-plane spin den- 
sity stated already below Eq. ([I]). For the spin density 
we find 



drs x 
d T s y 



Vp_ 

4 

Vp 



d y n 



d x n 



2t 
1 







o, 



(20) 
(21) 



4 ' 2r v 

from which we identify 2r as the spin relaxation time. In 



(20) and (21) we ignored the scalar electric potential, (j>. 

In order to compare with Ref. |16j we will now dis- 
cuss the diffusive limit. The particle density obeys the 
diffusion equation 



OtU - Ddj 



R' 



0, D= -^ Ttr ,T tr = 2T, (22) 



i.e., the particle current is j = —Dd R n. The transport 
time Ttr being twice as long as the scattering time r 
stems from the absence of backscattering in the helical 
metal. Equation (22) is consistent with Ref. jTHj, notice 
however the different definition of the diffusion constant. 
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Since the spin relaxation time is very short (it equals the 
transport scattering time) the spin dynamics is not diffu- 
sive. However, Eqs. (20 1 and (21) are still valid. In the 



diffusive limit the time derivative of the spin-density is 
small compared to the spin relaxation term, i.e. the spin 
density is given by the spatial derivative of the charge 
density so again we identify the general relation between 
particle current and spin density. 

We will now analyze the transport through the heli- 
cal metal when it is contacted via a tunnel junction to a 
ferromagnet. We will find an unconventional magnetore- 
sistance effect that arises since the tunneling probability 
between the helical metal and the ferromagnet is strongly 
angle dependent: the overlap of two spinors with polar- 
ization in rh- and rh'-directions depends on the angle 
between the two vectors, |(rh|rh')| 2 = | + • rh'. The 
states in the conduction band of the helical metal are po- 
larized as rh' = k x e z , so that the tunneling probability 
from a state rh in the ferromagnet into a state in the 
helical metal with momentum k depends on the angle 
between rh and k. 

In order to make these considerations more formal we 
introduce the tunneling Hamiltonian 



H t = 



d 2 xt(x)ijjl(x)il) F}S (x) +c.c, (23) 



where ipl (x) is the field operator for an electron with spin 
s in the helical metal and ipp a (x.) that for the ferromag- 
net. For a point-like (on the quasiclassical scale) contact 
at x = the tunneling amplitude is 



t(x) = ta(x). 



(24) 



Eq. (23) leads then to an additional contribution to the 



self-energy of the form 



E t (R,T;e) = |t| 2 5(R) 



d 2 k 



G F (R,T;k,e) (25) 



where Gf is the Green function of the ferromagnet. We 
assume that the conduction electrons in the ferromagnet 
can be described in terms of an incoherent superposition 
of majority and minority carriers, so we write 



Gf — GFtfPf + Gf,iPi, 



(26) 



where P-^j, = (ctq ± rh • er)/2 projects on states parallel 
antiparallel to rh. The kinetic equation ( 14 ) becomes 



(d T + vph ■ d R )go 



-go + - ( {go) + k • (if g ) 



- 7r|i| 2 <5(R)JV t (l + m . k x e z )(g Q - -(g t )) 

- vr|t| 2 5(R)^(l rh • k x e z )(g Q - *<^», (27) 

where g^^ and A^j, are the quasiclassical Green function 
and the density of states in the ferromagnet. The terms 
in the second and third line of this equation describe 




FIG. 1. Schematic view of system under consideration: elec- 
trons can tunnel between a ferromagnet (top) and a helical 
metal (bottom). The size of the contact is l x times l y , the 
current is assumed to flow only in x-direction. 



tunneling between the helical metal and the spin up or 
down band of the ferromagnet. As anticipated above, the 
tunneling probability between the helical metal and the 
ferromagnet is a k-dependent function. After the angular 
average and integrating over the energy we obtain the 
following continuity equation for the charge density, p — 
(-e)n, 



d T p + 9 R ■ j = -S(R) [G t (U - U t ) + G l (U~ UJ) 



- £(R)-^ — ie 2 x m-j, 
having introduced the tunneling conductances 



G 



7re' 



(28) 



(29) 



and the integrals 



U 



1 

2c" 



1 

4c" 



de(g n )- (30) 



The latter have the meaning of a voltage, or more pre- 
cisely of a magneto-electrochemical potential, cf. Eq. 
(18). Sources and sinks for the charge density in the he- 



lical metal appear on the right hand side of the equation 
due to tunneling. One observes that a tunnel current 
cannot only be driven by a voltage across the contact, 
but also by a current bias in the helical metal. 

For simplicity we will assume in the following equal 
magneto-electrochemical potentials for spin up and down 
electrons in the ferromagnet, t/^ = = Up Further- 
more, instead of a point-like tunnel contact we will con- 
sider an extented tunneling region as depicted in Fig. [T] 
This is achieved by replacing the 5-function on the right 
hand side of Eq. ( 28 ) by an appropriate function charac- 
terising the shape of the contact. In particular, we will 
study in detail a junction with a long extension l y i n y- 



directiorJ22l. We replace then the 5-function in Eq. ( |28[ ) 

by 



5(B) 



l~ 1 8(x) for - l y /2 < y < l y /2 
otherwise 



(31) 



4 



and we assume that current flows only in x-direction. 
The continuity equation thus becomes one-dimensional. 
Due to the remaining one-dimensional <5-function the cur- 
rent density jumps at x = 0. The size of this jump is 
determined by integrating the continuity equation with 
respect to the x-coordinate for a region close to the tun- 
neling contact with the result 



.]+ - 3- 



Gi j+ + j- : 



LN v F e 2 



(32) 

where j± = j(x = ±0) is the current right and left 
from the tunnel contact and we assumed that 5(x)j(x) = 
6(x)(j+ +j-)/2. When we fix the current left to the 
contact to zero we determine the tunnel resistance as 



Rt = 



1 



h Xprhy G^ — G 



T -I- 



G-f + G± 



2ly G^ + G 



(33) 



where for clarity we put back h. Through m y the tunnel 
resistance depends on the orientation of the ferromagnet 
with respect the the x-axis, similar to what was found 
by Burkov and Hawthorn^. The relative size of the 
magnetoresistance effect increases with increasing tun- 
nel conductance being controlled by the dimensionless 
parameters (G-f — Gi)/(e 2 /h) and Xf/I v - Apparently 
the tunnel resistance may even become negative for large 
enough values of these two parameters. However, as we 
will see below, when this happens equation (33 1 is no 



longer valid and a more careful treatment which takes 
into account a finite contact area of the tunnel junction is 
needed. To analyze a junction with a finite width l y and 
length l x we make the following replacement in Eq. (28): 
<5(R) — > (Ixly) when R is inside the tunnel junction 
and 5{YL) — > otherwise. This leads to the tunnel resis- 
tance 



Rt -G^TGl f 



Xf G^ - G_i 
l y e 2 /h 



with the function f(x) = x/(e x — 1). For small x, f(x) « 
1 — x/2 and one recovers Eq. (33). When x increases and 



becomes of order one, the resistance remains positive, but 
depends strongly on the sign of x, i.e., the orientation 
of the ferromagnet. For x > 0, the resistance decreases 
exponentially in the parameter {G^ — G^)/(e 2 /h){\F/l y ), 
while it increases linearly for x < 0. The tunnel junction 
is then acting as a spin-diode. 

In this final part of the paper we consider tunnel- 
ing from the helical metal into a normal, nonmagnetic 
metal. The tunnel resistance for the charge becomes 
R t = 1/(G^ + G^) since spin up and down have both the 
same density of states, — — iV-j-j, and also identical 
tunnel conductance. Nevertheless tunneling into a nor- 
mal metal is of interest since, as wc will demonstrate, the 
helical metal injects spin into the normal metal. Going 
through the same steps as in Eqs. ( 25 )— ( 28 ) but now for 



the spin-density in the normal metal we arrive at 
d tS a + d R -r = -^-^5(R)s a 



Rt N n 

h/2e 2 <5(R) 
Rt cNqvf 



(j x e,) OJ (35) 



where j a is the spin current in the normal metal and 
j the charge current in the helical metal. Whereas the 
tunneling of charge is controlled by the voltage across the 
tunnel junction, the tunneling of spin is controlled by the 
current density in the helical metal, i.e. a current flowing 
parallel to the junction. Thus injection of a pure spin 
current into the normal metal is possible. Furthermore, if 
spin tunneling in and out are balanced so that there is no 
spin current injection, there is a steady state determined 
by 



1 1 

N H eN a v F 



J x e 2 



which is equivalent to 



2N- 



n. 



normal metal 



N 



helical metal 



(36) 



(37) 



A normal metal on top of the helical metal thus amplifies 
the current-induced spin polarization. 

We finally turn to the question of how robust the ef- 
fects found are, having in mind that the effective Hamil- 
tonian ([T]) is only valid in the vicinity of the Dirac point, 
whereas the interpr etatio n of the experimentally ob- 
served Fermi surfaces ^ 5 ! 7 ! requires at finite doping terms 
that are quadratic and even cubic in kP^I. For example, to 
third order in k, the effective Hamiltonian for the surface 
states of the topological insulators Bi2Te3 and Bi2Se2 is^ 



H = wk x e z • er 



£ + + (38) 



(34) where the Dirac velocity contains a second order correc- 
tion, v — Wf(1 + oik 2 ), k± — k x ± iky, and a, m*,A are 
the parameters characterizing the strength of the higher 
order corrections. In order to understand the effect these 
extra terms in the Hamiltonian have on the magnetoresis- 
tance and the spin injection, recall that the origin of both 
is the helicity of the conduction electrons. The quadratic 
term in the Hamiltonian ( 38 1 has no spin structure, does 
not affect the helicity of the eigenstates, and therefore 
cannot qualitatively change our results. The cubic term 



on the other hand disturbs the helicity of the eigenstates 
(the angle between the velocity and the spin now de- 
pends on the position on the Fermi surface), so for a 
strong cubic term we leave the region where our results 
are reliable. 

In conclusion we have derived a kinetic equation for a 
helical model and have shown that the spin dynamics is 
constrained to follow that of the charge. In the diffusive 
regime for charge, the spin density is described in terms 
of the charge density gradient. Furthermore, when the 



5 



helical metal is placed in contact with a ferromagnetic 
metal, the tunneling current depends on the relative ori- 
entation between the current and the polarization in the 
ferromagnet. In the limit of large tunneling conductance, 
the device acts as a spin-diode. A helical metal in contact 
with a normal metal injects spin into the latter, with a 



rate that can be controlled independently from the injec- 
tion of charge carriers. 
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